Determination of subcatchment and watershed boundaries in a complex and highly urbanized landscape

Urban development significantly alters the landscape by introducing widespread impervious surfaces, which quickly convey surface run‐off to streams via stormwater sewer networks, resulting in “flashy” hydrological responses. Here, we present the inadequacies of using raster‐based digital elevation models and flow‐direction algorithms to delineate large and highly urbanized watersheds and propose an alternative approach that accounts for the influence of anthropogenically modified land cover. We use a semi‐automated approach that incorporates conventional drainage networks into overland flow paths and define the maximal run‐off contributing area. In this approach, stormwater pipes are clustered according to their slope attributes, which define flow direction. Land areas drained by each cluster and contributing (or exporting) flow to a topographically delineated catchment were determined. These land masses were subsequently added or removed from the catchment, modifying both the shape and the size. Our results in a highly urbanized Toronto, Canada, area watershed indicate a moderate net increase in the directly connected watershed area by 3% relative to a topographically forced method; however, differences across three smaller scale subcatchments are greater. Compared to topographic delineation, the directly connected watershed areas of both the upper and middle subcatchments decrease by 5% and 8%, respectively, whereas the lower subcatchment area increases by 15%. This is directly related to subsurface storm sewer pipes that cross topographic boundaries. When directly connected subcatchment area is plotted against total streamflow and flashiness indices using this method, the coefficients of variation are greater (0.93 to 0.97) compared to the use of digital elevation model‐derived subcatchment areas (0.78 to 0.85). The accurate identification of watershed and subcatchment boundaries should incorporate ancillary data such as stormwater sewer networks and retention basin drainage areas to reduce water budget errors in urban systems.     

Genetic variability of the Lessepsian migrant mussel Brachidontes pharaonis (Bivalvia: Mytilidae) in Tunisia

The present study used two mitochondrial markers (16S rRNA and COI) to assess the genetic diversity of a newly founded Lessepsian migrant mussel, Brachidontes pharaonis, in Tunisian waters. The species appears to be restricted to only one population in Rades Harbour, in the northern part of the country. Phylogenetic analyses revealed the monophyly of B. pharaonis in Tunisia. Both molecular markers revealed high genetic variability of the B. pharaonis population. Haplotype networks and demographic analyses confirmed the recent expansion events within this population. Multiple human-mediated introduction events involving several founder populations and intensive population growth rates are probably the main causes of the high polymorphism observed within this invasive mollusc.     

Nonlinear model to study filamentation instability of circularly polarized dispersive Alfvén waves in solar wind plasmas

This paper examines the small-scale solar wind turbulence driven in view of the Alfvén waves subjected to ponderomotive nonlinearity. Filamentation instability is known to take place for the case of dispersive Alfvén wave (DAW) propagating parallel to the ambient magnetic field. The ponderomotive force associated with DAW is responsible for wave localization and these webs of filaments become more intense and irregular as one proceeds along the spatial domain. The ponderomotive force associated with pump changes with pump parameters giving rise to different evolution patterns. This paper studies in detail the nonlinear evolution of filamentation instability introduced by dispersive Alfven waves (DAWs) which becomes dispersive on account of the finite frequency of DAW i.e., pump frequency is comparable to the ion cyclotron frequency. We have explicitly obtained the perturbation dynamics and then examined the impact of pump magnitude on the driven magnetic turbulence using numerical simulation. The results show steepening at small scales with increasing pump amplitude. The compressibility associated with acoustic fluctuations may explain the variation in spectral scaling of solar wind turbulence as observed by Alexandrova et al. (Astrophys. J. 674:1157, 2008).     

The spawning and early life-history of a New Zealand endemic amphidromous eleotrid,bluegill bully (Gobiomorphus hubbsi)

The larval stage of fishes is critical in determining their dispersal, survival and recruitment, but little is known of the larval behaviours and tolerances of amphidromous fishes, particularly in New Zealand. We report the results of a series of observational and experimental studies on bluegill bully (Gobiomorphus hubbsi), including spawning sites and behaviours, larval characteristics at hatch, phototactic responses of larvae, and larval survival at different salinity levels. Spawning primarily occurred in the lower reaches of the river, and larvae from different nests exhibited marked differences in, and trade-offs between, larval characteristics at hatch, potentially affecting larval success. Larvae were positively phototactic to intense light, an unexpected result based on diel drift patterns and international research. Finally, larvae exhibited markedly higher survival rates when reared at intermediate salinities compared to freshwater or seawater, suggesting estuaries may play an important role as nursery grounds for bluegill bully and other amphidromous fish.     

Otolith microchemistry indicates regional phylopatry in the larval phase of an amphidromous fish (Gobiomorphus hubbsi)

Knowledge of stock structure is key for the effective management of any fish species. Amphidromous fish, which live and spawn in freshwater but spend a pelagic larval period at sea, have typically been assumed to disperse widely during their larval phase, resulting in populations being sourced from a single unstructured larval pool. We used otolith microchemical analysis to examine the stock structure of bluegill bully (Gobiomorphus hubbsi), a declining amphidromous eleotrid endemic to New Zealand, along the west coast of South Island, New Zealand. Some drainages – even those in close proximity (c. 20?km) – were readily distinguishable based on otolith trace element concentrations, while little structure was evident between other geographically disparate locations. These results indicate that, at least in some cases, locally retained larvae, rather than a single unstructured larval pool, dominates recruitment. Management of bluegill bully and other amphidromous species must therefore consider the possibility of regionally distinct populations.     

On the coplanar eccentric non-restricted co-orbital dynamics

We study the phase space of eccentric coplanar co-orbitals in the non-restricted case. Departing from the quasi-circular case, we describe the evolution of the phase space as the eccentricities increase. We find that over a given value of the eccentricity, around 0.5 for equal mass co-orbitals, important topological changes occur in the phase space. These changes lead to the emergence of new co-orbital configurations and open a continuous path between the previously distinct trojan domains near the \(L_4\) and \(L_5\) eccentric Lagrangian equilibria. These topological changes are shown to be linked with the reconnection of families of quasi-periodic orbits of non-maximal dimension.     

Determination of Non-gravitational Effects in the Motion of Near-Sun Objects 321P, 322P, 323P,and 342P

At the beginning of this century, the SOHO space observatory discovered near-Sun comets with perihelion distances q ≈ 0.05 AU, which remained observable over several close encounters with the Sun. This became one of the surprises in studying the small bodies of the Solar System. Currently, there are objects that have already been observed in four (342P) and five (321P, 322P, and 323P) apparitions. In the present work, the estimates of nongravitational effects are obtained for these objects based on the pair-wise linkage of the apparitions. The calculations show that the observations of these objects are poorly represented if solely the gravitational forces are considered. The magnitude of nongravitational effects in the semimajor axis noticeably changes with time. The motion of all comets is significantly affected by the components of nongravitational forces that are perpendicular to the orbital plane.     

Wetlands and low-gradient topography are associated with longer hydrologic transit times in Precambrian Shield headwater catchments

The estimation of hydrologic transit times in a catchment provides insights into the integrated effects of water storage, mixing dynamics, and runoff generation processes. There has been limited effort to estimate transit times in southern boreal Precambrian Shield landscapes, which are characteristically heterogeneous with surface cover including till, thin soils, bedrock outcrops, and depressional wetland features that play contrasting hydrologic roles. This study presents approximately 3.5 years of precipitation and streamflow water isotope data and estimates mean transit times (MTTs) and the young water fraction (p_y) across six small catchments in the Muskoka-Haliburton region of south-central Ontario. The main objectives were to define a typical range of MTTs for headwater catchments in this region and to identify landscape variables that best explain differences in MTTs/p_y using airborne light detection and ranging and digital terrain analysis. Of the transit time distributions, the two parallel linear reservoir and gamma distributions best describe the hydrology of these catchments, particularly because of their ability to capture more extreme changes related to events such as snowmelt. The estimated MTTs, regardless of the modelling approach or distribution used, are positively associated with the percent wetland area and negatively with mean slope in the catchments. In this landscape, low-gradient features such as wetlands increase catchment scale water storage when antecedent conditions are dryer and decrease transit times when there is a moisture surplus, which plausibly explains the increases in MTTs and mean annual runoff from catchments with significant coverage of these landscape features.     

Changes in hillslope hydrology in a perched,shallow soil system due to clearcutting and residual biomass removal

Sustainable fuels legislation and volatility in energy prices have put additional pressures on the forestry sector to intensify the harvesting of biomass for “advanced biofuel” production. To better understand how residual biomass removal after harvest affects forest hydrology in relatively low slope terrain, a Before-After-Control-Impact (BACI) study was conducted in the USDA Forest Service's Marcell Experimental Forest, Minnesota, USA. Hydrological measurements were made from 2010–2013 on a forested hillslope that was divided into three treatment blocks, where one block was harvested and residual biomass removed (Biomass Removed), the second was harvested and residual biomass left (Biomass Left), and the last block was left as an Unharvested Control. The pre-harvest period (2 years) was 2010–11 and post-harvest (2 years) was 2012–13. Water table elevation at the upslope and downslope position, subsurface runoff, and soil moisture were measured between May–November. Mixed effect statistical models were used to compare both the before-after and “control” treatment ratios (ratios between harvested hillslopes and the Unharvested Control hillslope). Subsurface runoff significantly increased (p < .05) at both harvested hillslopes but to a greater degree on the Biomass Left hillslope. Greater subsurface runoff volumes at both harvested hillslopes were driven by substantial increases during fall, with additional significant increases during summer on the Biomass Left hillslope. The hydrological connectivity, inferred from event runoff ratios, increased due to harvesting at both hillslopes but only significantly on the Biomass Left hillslope. The winter harvest minimized soil disturbance, resulting in no change to the effective hydraulic conductivity distribution with depth. Thus, the observed hydrological changes were driven by increased effective precipitation and decreased evapotranspiration, increasing the duration that both harvested hillslopes were hydrologically active. The harvesting of residual biomass appears to lessen hydrological connectivity relative to leaving residual biomass on the hillslope, potentially decreasing downstream hydrological impacts of similar forestry operations.     

Solution for a plane strain rough‐walled hydraulic fracture driven by turbulent fluid through impermeable rock

The impact of turbulent flow on plane strain fluid‐driven crack propagation is an important but still poorly understood consideration in hydraulic fracture modeling. The changes that hydraulic fracturing has experienced over the past decade, especially in the area of fracturing fluids, have played a major role in the transition of the typical fluid regime from laminar to turbulent flow. Motivated by the increasing preponderance of high‐rate, water‐driven hydraulic fractures with high Reynolds number, we present a semianalytical solution for the propagation of a plane strain hydraulic fracture driven by a turbulent fluid in an impermeable formation. The formulation uses a power law relationship between the Darcy‐Weisbach friction factor and the scale of the fracture roughness, where one specific manifestation of this generalized friction factor is the classical Gauckler‐Manning‐Strickler approximation for turbulent flow in a rough‐walled channel. Conservation of mass, elasticity, and crack propagation are also solved simultaneously. We obtain a semianalytical solution using an orthogonal polynomial series. An approximate closed‐form solution is enabled by a choice of orthogonal polynomials embedding the near‐tip asymptotic behavior and thus giving very rapid convergence; a precise solution is obtained with 2 terms of the series. By comparison with numerical simulations, we show that the transition region between the laminar and turbulent regimes can be relatively small so that full solutions can often be well approximated by either a fully laminar or fully turbulent solution.